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ABSTRACT 



Context. Galaxy mass and environment are known to play a key role in galaxy evolution: studying galaxy colors as a function of 
redshift, galaxy mass, and environment offers a powerful diagnosis to disentangle the role of each. 

Aims. We study the simultaneous dependence of the fraction of blue galaxies f\,] ue on secular evolution, environment, and galaxy mass 
with a well-controlled cluster sample. We are thus able to study the evolution and respective role of the cessation of star formation 
history (SFH) in clusters caused by galaxy mass ("mass quenching") or by environment ("environmental quenching"). 
Methods. We defined an homogenous X-ray selected cluster sample (25 clusters with < z < 1 and one cluster at z ~ 2.2), 
having similar masses and well-defined sizes. Using multicolor photometry and a large spectroscopic sample to calibrate photometric 
redshifts, we carefully estimated f Uue for each cluster at different galaxy mass and cluster-centric distance bins. We then fitted the 
dependence of fti ue on redshift (z), environment (r/rjoo) and galaxy mass (M) with a simple model. 
Results. fbh,e increases with cluster-centric distance with a slope a = 1.2^2 3, decreases with galaxy mass with a slope y 
and increases with redshift with a slope y = 3.2*2,. The data also require for the first time a differential evolution with galaxy mass 
of fbiue with redshift, with lower mass galaxies evolving slower by a factor f = -4.1*H. 

Conclusions. Our study shows that the processes responsible for the cessation of star formation in clusters are effective at all epochs 
(z ^ 2.2), and more effective in denser environments and for more massive galaxies. We found that the mass and environmental 
quenchings are separable, that environmental quenching does not change with epoch, and that mass quenching is a dynamical process, 
i.e. its evolutionary rate is mass-dependent. Our study extends the downsizing-like scenario, where the most massive galaxies have 
their properties set at a very high redshift, to the cluster environment and all galaxies. It illustrates the need to disentangle galaxy mass 
and cluster-centric distance to properly estimate the behavior of f blue in clusters. 

Key words. Galaxies: clusters: general - Galaxies: clusters: individual: JKCS 041 - Galaxies: evolution - Galaxies: star formation 
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1. Introduction 

It is well-established that galaxy mass and e nvironment are key 
parameters in shaping galaxy properties (e.g., Butcher & Oemlen 
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Il978t iDresslerl Il980t iDe Propris et all I2003I) . At least in the 
local Universe, both change galaxy properties in a simi- 
lar d irection, for instance making them preferentially red der 
(e.g JVisvanathan & Sandagdl977r.lButcher & OemleJl978l) or 
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of ea rly-type morphology (e.g., iBlanton et al. I I2003I: iDressleil 
119801) . Furthermore, there is a correlation between galaxy mass 
and environment, because denser environ ments tend to be 
inhabited by mor e massive galaxies (e.g., iHogg et all 120031 : 



iBaldrv et atll20p6|). In the local Universe, the Sloan Digital Sky 
Survey (SPSS. lYork et al.ll2000h provided data from which one 
coul d quantify the role of each to an unprec edented precision 
(e.g jKauffmann et alJl2004tlPeng et a"Dl2010l) . 

However, when considering higher redshifts (z ~ 1-2), the 
analysis is complicated - and additionally there are less data 
available - by the disentangling of the precise role of secular 
galaxy evolution, of galaxy mass, and of environment on galaxy 
properties. Therefore, any analysis of galaxy evolution at z ~ 1- 
2 needs to simultaneously control for galaxy secular evolution, 
galaxy mass, and environment. Indeed, any bias in one of those 
three terms may confuse the interpretation of the obtained re- 
sult. For instance, a blueing related to a younger age at higher 
redshift may be mistaken for a blueing induced by a sampling 



biased preferentially toward lower density environments, which 
are richer in blue galaxies. 

Controlling for galaxy mass requires working with galaxy 
mass-selected samples, which presents the advantage of a bet- 
ter control of the sample. Indeed, galaxy mass has on average 
a much more regular evolution (increasing with decreasing red- 
shift) than galaxy rest-frame optical luminosity, which can tem- 
porarily increase with starbust, and then decrease. For example, 
a sample selected on rest-frame optical luminosity is very likely 
to be biased toward those temporary low-mass starbust galaxies 
at high redshifts, thus introducing a spurious increase in the frac- 
tion of blu e galaxies that disappea rs if a mass-selected sample is 
used (e.g jDe Propris et~a l. 2003). A more subtle question is the 
mass evolution of star-forming galaxies that will significantly in- 
crease their stellar mass with time. For example, an average star 
formation rate of 2 Mq yr~' (resp. 5 Mq yr~') for 5 Gyr (resp. 2 
Gyr) is sufficient to build 10 ll) M©. Therefore, selecting galaxies 
by stellar mass measured at the galaxy redshift may introduce 
a net inflow with decreasing redshift of galaxies that rise above 
the mass threshold adopted for the study. This uncontrolled in- 
flow with decreasing redshift of galaxies in the selected sample 
leaves the ambiguity between a real evolution of massive blue 
galaxies and a selection effect. To remove this ambiguity, one 
needs to control for mass not at the redshift of observation, but 
evolved at z — to select at high redshift the likely ancestors of 
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present-day galaxies of a g iven mass (e.g., lAndreon et al.ll2008t 
IRaichoor & Andreonll2012l) . 

Controlling for secular evolution requires acknowledging 
that the galaxy rest-frame color evolves with time. Indeed, galax- 
ies observed at higher redshifts will on average have bluer rest- 
frame colors than their local counterparts, because of their mean 
younger age and of the higher me an star formation activity of the 
Unive rse at higher redshifts (e.g.. lLiliv et al.llT996l iMadau et al.l 
fl998h . Using a non-evolving rest-frame color will classify an 
increasing number (with increasing redshift) of red-sequence 
galaxies as blue galaxies, thus possibly introducing a spuri- 
ous evolutionary trend. One can take this into account by us- 
ing evolving rest-frame colors to chara cterize galaxy proper- 
ties at different redshifts, as proposed by lAndreon et al.l (2006) 
and now adopted in several works (e.g., Haines et al.l 12009 ; 
iPeng et alJl2010l: IRaichoor & Andreonl2012l) . 



by analyzing a cluster sample spanning a wide redshift baseline. 
An observationally economic but still physically efficient proxy 
for measuring the star formation activity is galaxy color (e.g 



Controlling for environment raises the question of how to 
measure environment. Available data often oblige us to use a 
proxy for an environment estimate based on samples with either 
photometric redshift or that are rest-frame B-band luminosity- 
selected. However, at z ~ 1-2, these proxies may be either prone 
to large uncertainties or potentially biased toward (temporarily) 
blue overdensities. Again, this may introduce spurious trends 
or smooth existing ones (for instance, see discussion in §3 of 
lOuadri et al.ll2012l) . In this regard, clusters of galaxies are ideal 
laboratories. On the one hand, their cores unambiguously rep- 
resent the densest environment in the Universe at each epoch, 
where environmental processes consequently are the most effec- 
tive. On the other hand, they provide, at fixed redshift, handy 
samples of galaxies observed in similar conditions, with lit- 
tle contamination. Nevertheless, when using cluster-centric dis- 
tance as the environment's measurement in clusters, one should 
scale it with the cluster size (e.g. r2oo): if not, one risks to assign 
the same environmental measurement to different physical en- 
vironments, thus mixing age-related trend (redshift dependence) 
with environmental trend. The effect is of paramount importance 
because gradients in cluster populations are usually important. 
For instance, if one probes a cluster region within a fixed ra- 
dius (e.g. 1 Mpc), one will probe only central regions of large 
clusters but outer regions for small clusters, which possibly in- 
trod uces an artificial depende nce with respect to the cluster size 
(e.g. iMargoniner et al.ll2001l) . Moreover, even if challenging at 
z ~ 1-2, a robust cluster size estimation is needed to prevent 
any smoothing of pos sible environme ntal trends when co-adding 
data (see for instance Loh et al. 2008, and the use of the approx- 
imate cluster size estimate Bqc)- 

Additionally, if one uses galaxy clusters to study galaxy evo- 
lution, the cluster selection function should be as independent as 
possible of the (studied) galaxy population. As an example, by 
selecting preferentially clusters rich in blue galaxies at higher 
redshift (be cause they are easier to detect, as in the early 80's, 
lKronlll995h or richer in red galaxies ( because cluster s are de- 
tected thanks to the red sequence, as in lLoh et al.ll2008l with the 
RCS) likely introduces a selection effect whe n estimating the 
fraction of blue galaxies (lAndreon et al.ll2006h . It is preferable 
to select clusters by their X-ray emission, because at a given 
X-ray emission (Lx or Tx) clusters rich in blue galaxies are 
not fa vored/disfavored at a given cluster mass (lAndreon & Ettoril 
[1991 . Furthermore, it is well-known tha t cluster galaxy proper- 
ties depend on the cluster richness (e.g., Oemlerlll974l iDresslerl 
I1980L Dressi er et al 1 11 9971: lHansen et al.ll2009l) : a way to control 
this dependence is to select clusters with similar masses. 

In this paper, we aim at disentangling the role of galaxy secu- 
lar evolution, galaxy mass, and environment on galaxy evolution 



iButcher & Oemlerll984UWilliams et al.l2009HPeng et al.l201( 
We classified galaxies into two broad categories (blue or red) and 
defined the fraction of blue galaxies fbiue at a given galaxy mass 
M and scaled clu ster-centric distance rj r?no , in the vein of the pi- 
oneering work of IButcher & Oemlerl (11984k In the last decades, 
numerous works studied fpiue in clusters (e . g., De Propris et al.l 
120031: lAndreon et al.ll2006l lLoh et aUl2008t iHaines et al.ll2009l 
to name a few) but none of them studied the dependence of fu ue 
on redshift, galaxy mass, and cluster-centric distance at the same 
time. Analyzing simultaneously the dependence of fi,i ue on those 
three parameters allows one to put strong constraints on the var- 
ious processes that are responsible for the cessation of star for- 
mation activity ("quenching"). 

To achieve this, we considered the aforementioned remarks 
to our best ability. We used a sizable X-ray selected cluster sam- 
ple (25 clusters with < z < 1 and the JKCS 041 cluster at z ~ 
2.2) with well-defined r2oo and similar cluster masses, thus min- 
imizing the correlation between fbi ue and our cluster selection 
function. We then estimated fbiue as a function of galaxy mass 
M and scaled cluster-centric distance r/rjoo, defining blue/red 
galaxies considering also the stellar evolution of galaxies with 
time. To lead the analysis, we used a reasonable evolutionary 
model as a reference (exponentially declining star formation his- 
tory, fixed formation redshift), and then measured any deviation 
from this reference. This reference model does not aim to de- 
scribe the individual behavior of each galaxy, but to cancel out 
overall trends such as stellar aging or stellar mass increasing 
with time. Deviations from this reference model highlight the 
evolution and respective role of the quenching in clusters caused 
by galaxy mass or to environment. 

The plan of this paper is as follows. We define in Sect.|2]the 
cluster sample on which this study relies. Sect. [3] presents the 
data and their analysis, including our estimation of fbi ue - Our re- 
sults are presented in Sect.|4]and are summarized and discussed 
in Sect.gJ We adopt H Q = 70 km s _1 Mpc 1 , O m = 0.30, and 
Ha = 0.70 throughout. All magnitudes are in the AB system , 
corrected for Galactic extin ction usin g LSchlegel et al.l (119981) . 
Masses are computed with a IChabrierl (120031) initial mass func- 
tion and, if not stated otherwise, are defined by the mass of the 
gas that will eventually be turned into stars, i.e. corresponding to 
the integral of the star formation rate. 

2. Cluster sample 

We describe in this section our cluster sample. It was assem- 
bled from an X-ray selection and clusters with similar masses at 
different redshifts, using the cluster X-ray temperature as mass 
proxy. The X-ray selection ensures that our cluster sample selec- 
tion is unbiased by fbi ue , because the probability of inclusion of 
a cluster in the sample is independent of fbiue at a given cluster 
mass. At intermediate redshift, we relied on the X-ray selected 
cluster sample from the first five deg 2 of the X MM Large-Scale 
Structure survey (XMMLSS. lPierre et al.l2004l) . At low redshift, 
we extracted a sample from t he Highest X-ray FLUx Ga laxy 
Cluster Sample (HIFLUGCS, iReiprich & Bohringed 120021 an 
X-ray selected cluster sample, by selecting clusters with compa- 
rable cluster masses. We added to this sample the JKCS 041 clus- 
ter (z ~ 2.2. lAndreon et al.ll2009t lAndreon & Huertas-Companvl 
1201 lh . Though this cluster is different from the other ob- 
jects, because it is not X-ray selected, it is the highest red- 
shift cluster with a significant X-ray emission; moreover, its 
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of galaxies in clusters in the last 10 Gyr 
Table 1. Cluster sample 
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Fig. 1. Tx as a function of redshift for our cluster sample. 

mass/temperature is comparable to the rest of our sample. Figure 
Q]displays the temperature Tx of our cluster sample as a function 
of redshift. We note that our clusters (except for JKCS 041) have 
kT ~ 3 keV, typical of intermediate mass clusters. We have r2oo 
for all the objects, which is of paramount importance, because 
it is known that fi,i ue depends on cluster-centric distance. It has 
been derived from Tx or cr v as described in Appendix [A] The 
properties of the clusters are listed in Table [TJ 

2.1. XMMLSS sample (0.14 < z < 1.05) 

The main body of our cluster sample comes from the XMMLSS. 
This survey is located in the Wl area of the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS), which provides 
u*g'r'i'z' -band data. Mor e precisely, we selected all clusters 
from the first five deg 2 dValtchanov et al.1 12004 IWillis et alj 
12005b iPacaud et al]l2007l) . which are in the CFHTLS Wl field 
(see Figure [2j. We removed clusters with lower-than-average 
data quality: XLSSC011 and XLSSC021 (their galaxies satu- 
rate CFHTLS Wl data because they are too nearby), XLSSC 004 
(no Tx or cr v measurements, so we cannot estimate r2ooX 
XLSSC 017 and XLSSC 020 (located where CFHTLS Wl data 
are corrupted), XLSSC 022 (blended with clusters at similar 
redshifts). We finally assembled 18 clusters, 15 of which hav- 
ing rop q estimated from Tx measurements with the XMM tele- 
scope (IPacaud et al.l 120071) . 3 (XLS SC 007, XLSSC 14, and 
XLSSC 016) from cr v measurements dWillis et al.ll2005l) . 

2.2. HIFLUGCS sample (0.02 <z< 0.05) 

We completed our cluster sample with a low-redshift cluster 
sample, with similar masses (Tx). We started from the X-ray 
selected and X-ray flux-limited HIFL UGCS, all of which have 
Chandra temperature measurements dHudson et al.l 120101) . Out 
of the 64 clusters from the HIFLUGCS sample, we selected 
those with z spe c ^ 0.02 (to avoid shredding), kTx(keV) < 3.8 (to 
match cluster masses) and those falling in the SDSS field, thus 
obtaining seven low-redshift clusters. Their spatial positions are 
displayed in Figure [3] 

3. Data and analysis 

To measure fbi ue , we adopted a conservative approach which 
may return large error bars, but yields unbiased values. For in- 



Cluster 


Zspec kTx CT r 


'200 




(keV) (km s" 1 ) 


(Mpc) 


SDSS - HIFLUGCS 

iii/im n mn 1 A1+U.U4 nni 



A1367 

MKW8 

A2634 

A2052 

A2063 

A2657 



0.022 
0.027 
0.031 
0.035 
0.035 
0.040 



3.58 + o 
3.00+jj 
3.19 + !j 
3.35 + ° 
3.77+^ 



3 52 +,) 
J - J -0.11 



1.31 
1.17 
1.21 
1.25 
1.34 
1.28 



CFHTLS Wl - XMMLSS 



XLSSC 041 
XLSSC 044 
XLSSC 025 
XLSSC 027 
XLSSC 008 
XLSSC 013 
XLSSC 040 
XLSSC 018 
XLSSC 016 
XLSSC 014 
XLSSC 006 
XLSSC 012 
XLSSC 049 
XLSSC 007 
XLSSC 001 
XLSSC 002 
XLSSC 029 
XLSSC 005 



0.14 
0.26 
0.26 
0.29 
0.30 
0.31 
0.32 
0.32 
0.33 
0.34 
0.43 
0.43 
0.49 
0.56 
0.61 
0.77 
1.05 
1.05 



1.3+»-i 
2 + ^ 

1 O+0.6 

2.o41 



4.8 +0 - 
2.0 + ll 



9 9+0.9 
z..z._ 05 

3.2 + °i 
2.8 + | 
4.1 + 
3.7 + 



703!** 
416^ 



190+178 



0.66 
0.61 
0.80 
0.98 
0.60 
0.51 
0.67 
0.78 
1.25 
0.73 
1.26 
0.73 
0.75 
0.50 
0.88 
0.74 
0.79 
0.74 



JKCS 041 



2.2 7.3 + 



0.76 



Notes. T x are taken f rom | Hudson et aD {2010) for the 
SDSS/HIFLUGCS subsample, IPacaud et alj 120071) fo r the CFHTLS 
Wl/XMMLS S subsampl e and lAndreon et all d201 l"h for JKCS 041. 
o~ v are from IWillis et all d2005l) . r 2 oo are estimated from T x , except 
for for clusters XLSSC 007, XLLSC014 and XLSSC 016, which are 
estimated from cr v (cf. Appendix |Al. 



stance, to take into account the background (resp. stars), we 
adopted a minimum removal, that statistically accounts for the 
remaining background (resp. stars). 

To this aim, we needed to a) identify and remove stars; b) 
correct the underestimate of photometric errors listed in the orig- 
inal catalog, if any; c) correct for a (minor) residual photomet- 
ric offset; d) estimate redshift probability distribution functions 
p(z)\ e) mask potential contamination from other clusters at sim- 
ilar redshift and remove background sources; f) make sure that 
the different subfields have consistent colors, to prevent any sys- 
tematic offset in the estimated photometric redshifts. 

We describe in Sect. 13.11 the data used for the CFHTLS 
Wl/XMMLSS subsample, along with their analysis. The pro- 
cedure for the SDSS/HIFLUGCS subsample is similar, and is 
su mmarized in Sect. | 3.2I The a nalysis of JKCS 041 is described 
in iRaichoor & Andreonl We describe our estimation of 

fbiue in Sect. [33] 

We emphasize that we adopted a consistent method for all 
three cluster subsamples. 

3.1. CFHTLS W1 
3.1.1. Data 

For the CFHTLS Wl/XMMLSS subsample, we used the 
CFHTLS-T0006 data release of merged source catalogs 
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Fig. 2. CFHTLS Wl cluster subsample: spatial location of the 18 clusters of ourCFHTLS Wl subsample (red crosses). The magenta 
dashed lines represents CFHTLS Wl fields. For each cluster, blue circles have radii of 0.5x, lx and 2 x rooo- Gray areas represent 
masked regions for fbi ue estimation and background estimation, due to image corruption (rectangles) or potential clusters at similar 
redshifts as those of our cluster subsample for each field (disks and sectors). 



(u'g'r'i'z' bands, catalogs are available at the CFHT Science 
Data Ar chive site Eb- The VIM QS VLT Deep Survey project 
(VVDS, iLe Fevre et al.l 120051) gives spectroscopic redshifts, 
Zspec, °f several thousands of objects in the same area. Restricting 
the sample to z' < 22.5 (which corresponds to the faintest zf- 
band magnitude we work with), valid photometry and a secured 
Zspec (flag=3,4), yields a spectroscopic sample of 1652 galaxies 
(Zspec S 1) and 332 stars (z spec = 0). 

3.1.2. Star removal 

Stars spectroscopically identified by the VVDS cover less than 
1 deg 2 . Therefore, we chose the following method to homoge- 
neously remove stars in our sample. We adopted a conserva- 
tive app roach for star remova l by removing objects qualified as 
stars by ICoupon et al.l d2009l) according to their color and size 
that also have a photometric redshift incompatible with that of 
the cluster. Indeed, when cross-matched with objects that have a 
V VDS spectrum , 3% o f galaxies are wrongly identified as stars 
in ICoupon et al.l d2009l) . Although this is a low percentage, we 
additionally reduced the fals e identification rate by re-including 
objects classified as stars in ICoupon et al.l d2009l) . but verified 

that § Zd * 3x ^ Zrf p(z)dz > 0.9, where z c i is the cluster redshift and 

1 |http : / /wwwl . c ade- ccda . hia- iha . nrc- enre . gc . ca/ c£ht/| 
ITQQQ6 . html] 



p(z) is the photometric redshift probability distribution function 
output by Eazy (see Sect. 13.1. St . 

Stars not identified as such in this phase (12% of our VVDS 
stars) were removed at a later stage, during the photometric red- 
shift selection and the background statistical subtraction phase. 

3.1 .3. Photometric error correction 

We corrected the SExtractor dBertin & Arnoutsll996l l flux errors 
underestimation caused by the correlation of adjacent pixels dur- 
ing image resampling (e.g. | Casertano et al.l2000tlAndreonl200l|: 
iRaichoor & Andreonl 1201 21) as in iRaichoor & Andreonl (1201 2|) . 
and found a factor of 1.5, in agreement with previous studies 
(e.g.lllbert et al.l2006HCoupon et al.l2009T: IRaichoor & Andreonl 
120121) . 

3.1.4. Photometric calibration 

To prevent any systematic offset in colors between fields, which 
can lead to systematic offsets in the p(z) functions, we matched 
the color-color diagrams of stars falling in different fields. We 
selected a reference field (022539-041200), and bright stars 
(/' < 21, r2 < 3 pixels, S /N > 20 in all bands, terapix.flag 
< 1). We adopted ax 2 minimization with 3<x-clipping of stellar 
sequences in color-color diagrams, holding fixed all parameters 
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Fig. 3. SDSS cluster subsample: spatial location of the seven clusters of our SDSS subsample (red crosses). Magenta/red solid lines 
represent our SDSS field used for background estimation for each cluster (see Sect. I3.2t . For each cluster, blue circles have radii 
of 0.5 x, lx and 2 x r2oo- Gray areas represent masked regions for fuue estimation and background estimation, due to large-scale 
structures or potential clusters at similar redshifts as those of our cluster subsample for each field. 



Table 2. Photometric correction 6 applied to CFHTLS Wl sub- 
fields (m„ ew = m oU i + 5) 



Field 


u* 


g' 


(' 


z' 


022929-041200 


-0.092 


0.004 


0.001 


0.006 


022929-050800 


0.015 


0.011 


0.040 


0.015 


022539-050800 


-0.100 


-0.013 


0.007 


-0.021 


022150-041200 


0.009 


0.024 


0.039 


0.037 



Notes. We adopted 022539-041200 field and r'-band magnitudes as a 
reference. 



Table 3. Systematic offsets m/u - m meas between measured and 
best-fit model magnitudes for photometric redshift estimation 



Survey 


u 


g 


r 


i 


z 


CFHTLS Wl 


-0.07 


0.05 


0.01 


-0.03 


0.01 


SDSS 


-0.01 


-0.01 


0.02 


0.02 


-0.05 



but the intercept, which was frozen at the best-fit value of the ref- 
erence field. This relative photometric calibration (see Table [2} 
ensures that one measures unbiased p{z) all across our CFHTLS 
Wl fields. 



3.1.5. Redshift probability distributions p(z) 

To remove galaxies that are very probably in front of or behind 
the cluster (cf. j||3.1.6t . we estimated the full photometric redshift 
probability distributio n function p{z) for each galaxy, using Eazy 
(IBrammer et al. 2008) with default settings and an r-band mag- 
nitude prior (see Appendix IEV Systematic offset s in the photo- 
metric calibration between d ata and models (e.g.. lBrodwin et all 
120061: IBrammer et al.ll2008l) were estima ted with the same ap- 
proach as in iRaichoor & Andreonl (|20 1 2b : the found offsets are 
listed i n Table [3] These small s hifts, comparable with previous 
work s (lllbert etal.ll2006l 120091: ICouoon et al.ll2009t iBarro et al.l 
1201 U IRaichoor & Andreonl2012l) . were applied only during the 
p(z) estimation. 



3.1.6. Background removal 

We need to account for galaxies along the cluster line of sight, 
which we generally call background. We again took a conserva- 
tive approach by only removing galaxies that are very probably 
in front of or behind the cluster, and by statistically accounting 
for the residual background in a second step. 

First, we performed a photometric redshift selection by re- 
moving galaxies that are at z < Z c i - 0.05 x (1 + z c i) or z > 
Zd + 0.05 x (1 + Zd) at > 99 % confidence. This selection has 
been calibrated on our spectroscopic VVDS sample and was ac- 
complished by keeping objects with 



: rf -0.05x(l+z c () 



p(z)dz < 0.99 



and 



p(z)dz < 0.99. 



(1) 



(2) 



c rf +0.05x(l+z rf ) 



Figure |4]illustrates such a selection criterion for three galaxies. 

This selection is conservative, because it keeps galaxies in 
the sample that have a low probability to belong to the cluster. 
Eqs.(Q]i and (0 wrongly remove < 1% of galaxies from our spec- 
troscopic VVDS sample (13/1652). To illustrate the efficiency of 
this background removal, we display in Figure |5]the cumulative 
p(z) function for all objects in the background area for four clus- 
ters, before and after applying the selection described by Eqs.fUJi 
& (0. We observe that this selection removes > 50% of back- 
ground objects. We note that the same amount of objects are 
removed for cluster outer regions (1 < r/r2oo < 2), but less for 
cluster core regions (r/r^oo < 1)> because the background con- 
tamination is proportionally lower towards the cluster center. 

We thus observe that Eqs.fUJi & (0 fulfill our conservative 
approach criterion: it significantly reduces the sample, with re- 
jecting a non-significant part (< 1%) of the galaxies we are in- 
terested in. Changing the {0.05, 0.99} values toward a more strin- 
gent selection could reduce the contamination in our sample, but 
at the cost of removing a more significant part of the galaxies we 
are interested in. Changing the {0.05,0.99) values toward a less 
stringent selection would only increase the contamination in our 
sample. 

As a second stage of the background subtraction, we aimed 
to select independent control fields taken from the very same 
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Fig. 4. Illustration of the selection criterion applied in Eqs.([T} 
and (f2]i with z c i = 0.61: we display p(z) for three galaxies. For 
each galaxy, we report above the value for I\ (resp. h), the inte- 
gral in Eq.(Q~|i (resp. Eq.©), and whether the galaxy is selected 
or not. Vertical green solid and dashed lines represent z c i and 
Zd ± 0.05 x (1 + Zd), respectively. 
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Fig. 5. Cumulative p(z) function for all objects in background 
area for four clusters, before (black curves) and after (ma- 
genta curves) applying the selection described by Eqs.([TJ & 
©. Vertical green solid and dashed lines represent z c i and 
Zd ± 0.05 x (1 + Zd), respectively. 



data and CFHTLS Wl field as the cluster galaxies, to guaran- 
tee homogeneity across the samples. Furthermore, because we 
aimed to stack data for different clusters, we also needed inde- 
pendent data for the background, so that we did not overstate the 
quality of the control sample. 

We selected for each CFHTLS Wl field a control sample in 
the same field in the following way. First, we removed objects 
within 2 X 7-200 around each cluster and within circles (0.05° ra- 
dius) around potential clusters at similar redshifts (gray areas in 
Figure 0. The several studies on clusters in the CFHTLS Wl 
fields ensure that contaminating clusters are taken into account. 
Then, we randomly divided the remaining objects into four sub- 
sample to reduce the effect of any potential structure, which as- 
sembled four independent control subsamples (each with a sur- 
face of ~0.2 deg 2 ). We randomly associated each cluster to a 



control subsample. When clusters in a common field had simi- 
lar redshifts (always less than four in our sample), we associated 
them with different control subsamples in this field, thus never 
attributing the same background twice. 

In addition, this step also subtracts any star that has not been 
identified as one based on its colors. 



3.2. SDSS 

SDSS data were analyzed identically. For the SDSS/HIFLUGCS 
subsa mple, we worked w ith data in ugriz bands from the SDSS 
DR8 dAihara et alj|2011l) . We used the ModelMag for galaxy col- 
ors, and cModelMag and PsfMag as an approximate total mag- 
nitude for galaxies and stars, respectively. 

Our SDSS galaxy sample is selected from the Galaxy view 
of the SDSS DR8, to which we applied a cut in magnitude 
(cModelMagjr < 19), a spatial cut (10° x 10° around each 
cluster) and a criterion to reject saturated stars misclassified as 
galaxies (FiberMag.u > 22 - 0.5x ModelMag_u). By matching 
this catalog with the spectroscopic informations from the main 
galaxy sample (petroMag_r < 17.8 in SpecPhoto view), we 
obtained ~30,000 spectroscopic redshifts (z S pec 5= 0.4), which 
we used to assess the quality of our background removal. 

To check that the photometry does not present any system- 
atic offset between our fields, we used more than 10,000 spec- 
troscopic stars from the SpecPhoto view located in our fields in 
in the same way as for the CFHTLS Wl fields (Sect. 13X41 We 
found indeed no offset greater than 0.01 mag for all filters and 
all fields. 

We estimated the p(z) in in the same way as for the CFHTLS 
Wl fields (Sect. [3~T31 l. We applied the photometric offsets listed 
in Table E] 

For background removal, we selected galaxies as in Sect. 
13.1.61 This criterion removes ~20% of the background objects, 
and wrongly removes ~0.05% of the spectroscopic objects. For 
each cluster, we selected as a control sample objects within a 
10° x 10° area around each cluster, excluding objects within 
2 x r2oo around each cluster and within circles (0.3° radius) 
around potential clusters at similar redshifts. For clusters A2052 
and A2063, we defined two independent control subsamples as 
described in Sect. 13.1.61 and using the region within the red solid 
line in Figure [3] The resulting control samples are ~90 deg 2 
(resp. -50 deg 2 ) for the MKW4, A1367, MKW8, A2634, and 
A2657 clusters (resp. the A2052 and A2063 clusters). 

3.3. Settings for the estimation of fu M 

We describe in this section our procedure for consistently mea- 
suring the fraction of blue galaxies fu U e- The procedure de- 
scribed here is similar to the one followed in Section 4 of 
iRaichoor & Andreonl d2012b for JKCS 04 1 . 

For each cluster, we considered three radial bins (annuli) at 
distinct cluster-centric radii, defined by r/r^oo ^ 0.5 , 0.5 < 
r / r 200 - 1, and 1 < r/r^oo - 2. 



3.3.1. Probed color 

To probe the same rest-frame color at all redshifts as consistently 
as possible, we worked in a color-magnitude diagram that clos- 
est corresponds to the (u - r) vs r rest-frame, thus probing the 
4000 A break. Figure |6]illustrates the filters chosen at each red- 
shift of our sample. Because our reddest band for the CFHTLS 
Wl sample is the z band, this implies a drift in the probed rest- 
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Fig. 6. (u - r)-like rest-frame probed color: for each cluster of 
our sample, we plot the rest-frame probed color as a function 
of the redshift. The shaded areas indicate the wavelength inter- 
vals where the filter response convolved with a typical elliptical 
galaxy spectrum is above half of its maximum value. The verti- 
cal bars represent the same for the redshifts of our clusters and 
the corresponding chosen filters. 



frame red band when the redshift is approaching 1 (z — 0.77 and 
z = 1.05): we are here limited by the data, but s till improve on 
previous works (see for instance Loh et al. 2008). 

The top panel of Figure ^illustrates another possible choice 
of filter pairs, matching the (B - V) rest-frame color. The bottom 
panel shows that f\,\ ue - estimated as described below - does not 
depend on the precise sampling of the filters, given that they 
bracket the 4000 A break. Our results are robust regarding the 
rest-frame color: using the (B - V) color instead of (u - r) leads 
to similar values of fbi ue (see Figure|7]i. 



3.3.2. Mass bin defintion 

Fo r our values of mass M w e refer, as in previous works, to 
fhe lBruzual & Chariot) (120031) model (2007 version, CB07 here- 
after) mass, and specifically to the mass of the gas that will 
eventually b e turned into stars, i. e . the integral of the star for - 
mation rate dAndreon et al.l 120081: iRaichoor & Andreonll2012l) . 
This mass definition has the advantage to compare galaxies at 
each epoch that will end up with the same stellar mass, which 
means that we include in each mass bin the same galaxies at high 
and low redshifts. In the color-magnitude diagrams, the shape of 
these loci of constant mass differ from that which corresponds to 
the mass in stars at the redshift of observation (M :o b s ) on the blue 
end (see Figure|8]i: the higher the redshift, the more fainter/bluer 
objects are included in our definition. For instance, a galaxy with 
an exponentially declining SFH oc exp[-(f/r)] with r = 3 Gyr 
would have a mass in stars of ~30/50/90% of its final mass in 
stars after ~0. 8/1. 6/5. 3 Gyr: had we chosen the mass in stars 
at redshift of observation as a definition of stellar mass, such a 
galaxy would not be included in our sample at high redshift, then 
would enter it at a given epoch and drift from one mass bin to 
another at an irregular pace, thus biasing the result. We discuss 
in Section l5T2l the impact of this definition on our results. 

We used four galaxy mass bins (log(M/M0) in [9.92, 10.33[, 
[10.33, 10.73[, [10.73, 11.13[, and [11.13, +oo[). We hereafter 
refer to those g alaxy mass bins by using the mean value of 
each bin (using a lSchechterll 1 9761 function), i.e. (log(M/Mo)) ~ 



p 20 




Fig. 7. Dependence of fbi„e on the rest-frame probed color. Top: 
like Figure|6] but for B — V rest-frame color. Bottom: we compare 
the estimated fi,[ ue when using the rest-frame B—V and u — r col- 
ors for 84 estimated values (10 clusters for which are available 
independent filter pairs probing the rest-frame B - V and u — r 
colors, 1 to 4 mass bins, 3 radial bins). We stress that some in- 
dividual estimates of fbi„ e have large error bars (see Figure HOl : 
the 68% shortest intervals for fbi„ e (B - V) and fbi ue {u - r) all 
intersect. 



10.14, 10.54, 10.94, and 11.47. These (model) masses are com- 
puted for solar metallicity, a formation redshift of z form = 5 (set- 
ting Zform to 4 or 6 does not change our results), and either SSP or 
an exponentially declining star-forming t model with < SFH t 
(Gyr) < 10. The mass values used to define mass bins correspond 
to the mass of an SSP model having at z — a V-band absolute 
rest-frame magnitude M v of -17.8,-18.8,-19.8, and -20.8. 

The depth of the data limits the sampled mass range: for each 
cluster, we required that the lowest mass cut is brighter than a 
signal-to-noise ratio of 5 in the considered color-magnitude dia- 
gram. This constraint - probing higher masses at higher redshifts 
- restricts the mass range probed, but ensures that we work on 
mass-complete samples and allows a direct comparison of fbi m 
values at different redshift and mass bins. 

3.3.3. Blue/red definition 

We define a galaxy as blu e if it is bluer than a CB07 mode l 
with r = 3.7 Gyr, as in lAndreon et ail (120041 120061 l2008l> . 
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Fig. 8. Illustration of a mass bin definition (z bs — 0.49, 10.73 < 
log(M/M0) < 11.13): thick magenta solid lines represent the 
loci of constant mass M used in this study (the mass evolved to 
z = 0), thick blue dashed lines the loci of constant mass M~„/„ 
in stars at z bs> and gray shaded areas represent the difference 
between those two definitions. The horizontal dotted line repre- 
sents the threshold used to define blue/red galaxies. 

lLoh et alj ( 120081) . and lRaichoor & AndreorJ(l2012l) . This galaxy 
will be bluer by 0.2 mag in B - V than red-sequence galaxies 
at z — (which would be a blu e galaxy by the original defini- 
tion of iButcher & Oemlerl [19841 A s a compar i son, o ur defini- 
tion agrees at z — with that used in lPeng et al.l (l2010h . and effi- 
ciently splits the red-sequence and the blue cloud of their SDSS 
data. The rationale behind our choice is to take into account the 
stellar evolution of galaxies with time. Figure [9] illustrates the 
rest-frame u - r color evolution for different exponentially de- 
clining r-SFHs: the different tracks do not cross each other. In 
other words, if a galaxy follows an exponentially r-SFH with 
t < 3.7 Gyr (resp. t > 3.7 Gyr) during its lifetime, its color ac- 
cording to our definition will remain red (resp. blue). However, 
we note that this figure assumes that all galaxies have the same 
Zform- We will return to this point in Sect. [5] 

In addition, we stress that the independence of fi,i ue on the 
probed rest-frame color (Sect. 13.3. It is a c onsequence of our 
choice to define red and blue galaxies (see iFairlev et al.ll2002l 
for an opposite result when using the red-sequence to split red 
and blue galaxies). 



4. Results 

We present in this section the dependence of fbi m with galaxy 
mass, cluster-centric distance and redshift. We recall that the 
evolution of fu ue with redshift traces the fraction of cluster 
galaxies that have stopped forming stars, under the assumption 
of an average exponentially declining star formation history for 
cluster galaxies. Indeed, under our working assumptions and our 
definition of blue/red galaxies, any quenching implies a lowering 
of fbiue with redshift, whereas no quenching implies an fbi ue con- 
stant with redshift. Hence our results present the evolution with 
redshift of the star formation activity in clusters as a function of 
galaxy mass and environment. 

Before performing a general analysis, we start by a simple 
analysis that makes partial use of the data (Sect. 14. Il l and that 
therefore does not show trends in a similar clear way (Sects. 14.21 



Fig. 9. Model rest-frame u — r color evolution for different expo- 
nentially declining t-SFHs (CB07 models, solar metallicity and 
Zform = 5). According to our definition a galaxy is classified as 
red (resp. blue) if it is redder (resp. bluer) than a model with 
t = 3.7 Gyr. The vertical arrows indicate the redshifts of our 
cluster sample. 

& I4.3I I. Our main result is presented in Section POI and illustrated 
in Figure [13] environmental and mass quenching are separable, 
environmental quenching does not change with epoch, and mass 
quenching is a dynamical process, i.e. its evolutionary rate is 
mass-dependent. 

4. 1 . Results for individual clusters 

For each cluster, galaxy mass bin, and radial area we com- 
puted the blue fraction fbi ue that accounts for residual back- 
ground galaxies (i.e. along the line of sight, kept by the criteria of 
Eqs.(Q]i & (0, and not belonging to the cluster) using our control 
samples (as defined in Sect. 1 3 . 1 . 6ft and foll owing the Bayesian 
methods introduced in lAndreon et al.l d2006). We adopt uniform 
priors for the parameters. 

Figure[10]displays the individual fi,i ue profiles for our cluster 
sample, grouped by redshift and mass bins. We defined our red- 
shift bins so that each bin spans a period of ~2 Gyr (0 < z < 0. 1 6, 
0.16 < z < 0.37, 0.37 < z < 0.65, 0.65 < z < 1.1, and z ~ 2.2). 
We hereafter refer to those redshift bins by using the median 
value of our cluster sample in each bin, i.e. (z) ~ 0.03, 0.31, 
0.49, 1.05, and 2.2. We can already observe some indications of 
the trends that we later perceive in stacking analysis. 

We observe a different behavior for different galaxy mass 
bins: at (z) ~ 0.31, fuuJj '/ '^oo ^ 0.5) decreases (from ~0.4 to 
~0) with increasing galaxy mass. fbi ue seems to increase with 
increasing r/r2oo: for instance, at (z) ~ 0.03 for the less massive 
galaxies, fbi ue increases with increasing cluster-centric distance. 
For the most massive galaxies, fbiuJj I ^200 < 0.5) ~ at all 
redshifts. 

However, the unavoidable noise - because of the finite and 
usually small number of member galaxies per cluster at a given 
mass and cluster-centric distance bins - makes it difficult to draw 
a general picture. 

4.2. Results for stacked clusters 

We now stacked our sample according to redshift and mass bins 
to strengthen the estimated fbiue- F° r each panel of Figure [TO] 
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Fig. 10. fbiue for individual clusters as a function of cluster-centric distance {r/rooa) for different bins of redshift (increasing right- 
ward) and galaxy mass (increasing downward). The computed value of fbiue already takes into account the mean aging of stars with 
increasing redshift by defining the color compared to an exponentially declining SFH model with t = 3.7 Gyr. Error bars represent 
the shortest interval including 68% of the possible posterior values for fbiue, and are plotted as a dashed line when this interval is 
larger than 0.66, indicating that fbiue is not constrained at all. Radial cluster-centric bins are indicated by gray filled triangles on the 
x-axis and horizontal gray dashed lines indicate the minimum and maximum allowed values for fbiue- 



we computed the mean fbi ue for all clusters belonging to this 
panel, multiplying the individual data likelihoods, and then us- 
ing the Bayes theorem as in the previous section. We recall that 
the combined data are independent, in particular the control sam- 
ples. The values of fbi ue are reported in Tableland displayed in 
Figure [TT] as dots with error bars. Now the trends are emerging 
quite clearly. 

Regarding the dependence of fbi ue with cluster-centric dis- 
tance, our data seem to indicate an increase in fbiue with increas- 
ing r/r2oo, for most redshift and mass bins. The cluster color 
profile at given galaxy mass tends to flatten with decreasing red- 
shift. In particular, the flat profile is achieved at the lowest red- 
shift and most massive galaxies. For all mass and radial bins, our 
data suggest a decrease of fbiue with decreasing redshift. For our 



two lowest mass bins, this evolutionary trend is strong. For the 
highest mass bin, this decrease is marginal, because fbiue already 
has very low values at high redshifts. 

Finally, we also see a trend when looking at the dependence 
of fbiue with galaxy mass. Indeed, for our two lowest redshift 
bins, where our data span four mass bins, we observe a clear 
trend at all radial bins: at fixed redshift and cluster-centric dis- 
tance, fbiue decreases with increasing galaxy mass. For instance, 
for the innermost radial bin and (z) ~ 0.03, fbiue decreases 
continuously (with increasing galaxy mass) from 0.11*221 to 
O-OI^qq}. For (z) ~ 0.49, though our data only span two mass 
bins, we observe the same trend. Our analysis shows that there 
is - in addition to the known mean aging of the stellar popula- 
tions - a decrease at all radii in fbiue with decreasing redshift, 
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Fig. 11. fbiue for stacked clusters as a function of cluster-centric distance (r/ r2oo) for different bins of redshift (increasing rightward) 
and galaxy mass (increasing downward). The computed value of fu U e already takes into account the mean aging of stars with 
increasing redshift by defining the color compared to an exponentially declining SFH model with t = 3.7 Gyr. Error bars represent 
the shortest interval including 68% of the possible posterior values for fbiue- Radial cluster-centric bins are indicated by gray filled 
triangles on the x-axis and horizontal gray dashed lines indicate the minimum and maximum allowed values for fbu, e - Yellow 
shaded areas with a solid blue line represent the posterior mean and 68% confidence interval of the modeling of Eq.fO, fitting the 
225 individual f\,\ ue measurements. Yellow slanting hatched areas with a dashed blue line represent the prediction/extrapolation of 
this model for bins where we do not have data. Cyan vertically hatched areas with solid magenta line (in two panels only) illustrate 
how the model fit fails if the £ term in Eq.[3]is not included. 



the intensity of which decreases when the galaxy mass increases, 
and disappears for the most massive galaxies. 



We note that in one panel of Figure QT| ((z) ~ 0.31 and 
<log(M/MQ)) ~ 1 1 .47), our data present an odd behavior for 
?7?"2oo ^ 1: the value of fu ue is quite high and does not reflect 
the general trend, as the mismatch with the fit (detailed below) 
illustrates. However, we keep in mind that the value is just at 

2cr. 



4.3. Results for the whole sample 

We now fitted the 225 values of fu U e (26 clusters, 3 radial bins, 
1 to 4 galaxy mass bins) at once by modeling at the same time 
the dependence of fbi ue on galaxy mass M, redshift z and cluster- 
centric distance r/r2oo- In this section, we are no longer ignoring 
redshift differences inside each panel, and the collective use of 
the whole dataset allows us to strengthen more the trends even. 

As detailed below, to described the data, we also need an 
interaction term between galaxy mass and redshift, a term that 
allows galaxies of different mass to evolve at different rates. We 
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Table 4. fi,i ue for our cluster sample stacked by redshift and galaxy mass bins. 
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Notes. The values of f blue are those plotted in Figure QTJ 



use the following model: 

fbiue(r/r 2 oo,M,z) = ilogit A + 

a ■ log (r/(0.25 ■ r 200 )) + 
/3 ■ (log(M/M G ) - 11) + 
y - (z - 0.3) + 

f .aog(M/Mo)-ll)-(z-0.3) 



(3) 



where ilogit(x) = (1 + exp(-x)) -1 ensures that < fu U e ^ 1. 
We adopt uniform priors for the parameters Aq, a, f3, y, and (. 

The model fit results are plotted in Figure QT| as solid blue 
lines and yellow shaded areas (68% confidence interval), while 
the model prediction/extrapolation is plotted as dashed blue lines 
and yellows hatched areas. 

The only motivation behind the chosen parametrization is the 
adoption of an additive model that fits our data and ensures that 
< fbiue - 1 ■ The adopted model is the simplest one fulfilling 
those criteria. If we do not include the f term, the model can- 
not fit data in the redshift bins (z) ~ 0.31 for galaxies within 
our two lowest mass bins (see cyan vertically hatched areas in 
Figure flTT i and, to a lesser extent, in the panel (z) ~ 0.49 and 
(log(M/M0)) ~ 10.94. In a similar manner, if we replace the 
term crossing M and z by a term crossing r and z or r and M, the 
model fails in fitting the data. Furthermore, we note that adding 
a term crossing r and z (in addition to our £ term) does not im- 
prove the quality of the fit. The f\,\ ue evolution -with redshift is 
therefore more sensitive to the galaxy mass than to its cluster- 
centric distance. 

The posterior probabilities for the model parameters are pi ot- 
ted in Figure [12] which shows how our knowledge about the 
parameters changes from before (prior, blue line) to after (poste- 
rior, black/red lines) the data. In short, the posterior distribution 
of all parameters is much more concentrated than the prior, i.e. 
data are highly informative about these parameters and conclu- 
sions on these parameters do not depend on the adopted prior. 

The values agree with our qualitative trends put forth in 
Sect. 14.21 with fbi„e increasing with cluster-centric distance 
(a = 1.2+q 4 ), decreasing with galaxy mass (J3 = -3.8+jjj), 
and increasing with redshift (y = 3.2+°^). ( measures the rel- 
ative speed - when compared to the evolution of galaxies with 



log(M/M0) = 1 1 - at which f\,\ ue evolves as a function of galaxy 
mass. Because £ is negative (f = -4.1+q '), galaxies with lower 
masses will on average evolve on longer timescales. We note that 
our model is constrained at z » 1 only by the JKCS 041 cluster 
(if we repeat the fit without the JKCS 041 data, the model ap- 
pears to be entirely unconstrained in the (log(M/MQ)} ~ 11.47 
and z ~ 2.2 window). 

To better visualize the dependence of fu U e on those three pa- 
rameters (M, r/>2oo, and z), we display in Figure [13] the vari- 
ations of fbiue on two of them when fixing the third one. The 
information (data and model fit/predictions) is exactly the same 
as in Figure [TT] only presented differently. We recall that our 
model is only constrained where there are some data (shaded 
areas), and that the fit was performed on the individual 225 val- 
ues of fbiue- This figure illustrates on the one hand that at fixed 
cluster-centric distance, ft,i„ e evolves at different paces for differ- 
ent galaxy stellar masses (panels a and b), captured by the ( term 
in our model. On the other hand, we can see in panel c that the 
model evolution of fbiue with z at fixed galaxy mass has a similar 
shape for our values of r/r2oo- Lastly, panel d summarizes the ra- 
dial profiles of fbiue at a fixed redshift for different galaxy mass 
bins. 



4.4. Results for the original Butcher & Oemler cut 

Adopting a mass thresh old corresponding to the original 
iButcher & Oemlerl (11984 one (log(M/M ) > 10.53), we ob- 
serve no evolution of fbiue with redshift for our data (z < 0.43) 
(see Figure HM. However, when extrapolating the model to z = 

1 .05, we observe an inc rease in fbi„e, in broad agreement with the 
measurement done by lAndreon et al. (2008| for the RzCS 052 
cluster. 

That almost no evolution is found at z < 0.43, whereas we 
do detect one for (log(M/M0)} ~ 10.54, nicely illustrates the 
benefit of splitting our sample into fine galaxy mass bins. 



5. Discussion and conclusion 

Our aim was to put constraints on the different processes respon- 
sible for the cessation of star formation activity ("quenching") in 
clusters through a careful study of the dependence of the fraction 
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Fig. 12. For each parameter fitted with the model of Eq.(|3), we plot as a black lines the posterior probability function and as yellow 
shaded area the shortest interval including 68% of the possible posterior values. We report the point estimate and the 68% confidence 
interval on each panel, and overplot as a thick red line the corresponding Gaussian curve. The blue lines represent the uniform prior 
used for the fit. 




Fig. 13. Dependence of fbi ue on M and z at fixed r/r2oo (panels a and b), on z and r/r2oo at fixed M (panel c), and on M and r/r2oo 
at fixed z (panel d). Data points with error bars are the same as in FigureQl] Shaded areas represent the posterior mean and 68% 
confidence interval of the modeling of Eq.(|3]i, fitting the 225 individual fw m measurements. Prediction/extrapolation of this model 
for bins where we do not have data are plotted as hatched areas and dashed lines. X-axis bins are indicated by gray filled triangles 
on the x-axis and horizontal gray dashed lines indicate the minimum and maximum allowed values for fbiue- 



of blue galaxies fbi„e in clusters that span a broad redshift base- 
line. We have built a cluster sample consisting of seven clusters 
(0.02 < Zspec < 0.04) from the HIFLUGCS survey, eighteen 
clusters (0.14 < z spec < 1-05) from the XMMLSS survey, and 
JKCS 041 (z ~ 2.2). Because the HIFLUGCS and XMMLSS 
surveys are X-ray selected, our cluster sample is unbiased re- 
garding the fraction of blue galaxies fbi„ e at a given cluster mass; 
in addition, our clusters were chosen to have similar masses. For 
different galaxy mass and cluster-centric distance bins, we care- 
fully estimated fbiue for each cluster and fitted the dependence of 
fbiue with redshift, galaxy mass and cluster-centric distance. We 
did not attempt to invert the noisy data because of the degeneracy 
between Zf or m, age (and metallicity), but instead we modified the 
model until it fitted the data. The reference model has an expo- 
nentially declinining SFH, and hence an evolving mass (if stars 
are formed, mass in stars cannot be fixed). Deviations from this 
reference model highlight the evolution and respective role of 
the quenching in clusters due to galaxy mass or to environment. 

Our main result is displayed in Figures QT| and [13] We re- 
call that in those figures we have already taken into account that 
galaxies are on average younger at higher redshifts, by using an 
exponentially declining T-model to define whether a galaxy is 
red or blue. This means that any measured evolution of f\,\ ue with 
redshift reflects a phenomenon in addition to the younger mean 
age of the Universe and the secular increase in the star formation 
rate. We also recall that our model fitting was performed on the 
225 individual measurements of fu m . Moreover, our definition 



of galaxy mass ensures that, for each mass bin, the selected high- 
redshift galaxies are those that would end up in the low-redshift 
galaxies selection. 

We are therefore able to disentangle for the first time the role 
of secular evolution, galaxy mass and cluster-centric distance 
on galaxy evolution in clusters. We found that fbiue decreases 
with decreasing redshift z, with decreasing cluster-centric dis- 
tance r/r2oo, and with increasing galaxy mass M. This means 
that the processes responsible for the cessation of star formation 
in clusters are effective at all epochs (z < 2.2) and more effec- 
tive in denser environments and for more massive galaxies. Our 
(simple) modeling of the data shows that the dependence of f,i ue 
on galaxy mass evolves with redshift, while this is not the case 
for the dependence of fbiue on cluster-centric distance: the inten- 
sity of mass quenching evolves with redshift at different paces 
for different galaxy masses, whereas the environmental quench- 
ing does not evolve with redshift. Additionally, our data did not 
need any correlation between galaxy mass and cluster-centric 
distance in the model fitting, which means that mass quenching 
and environmental quenching are separable. 

5. 1. Remark on taking into account secular evolution when 
estimating fbiue 

We illustrate in this section the advantage of taking into account 
secular evolution w hen estimating fbiue- using a non-evolving 
color threshold (e.g. iPresotto et al.ll2012h will result in using a 
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Fig. 14. Evolution of fbiue for individual (top panels) and stacked (bottom panels) clusters as a function of cluster-centric distance 
(r/r? M)) for different bins of red shift (increasing rightward) for a mass-selected sample with a mass corresponding to the original 
My in lButcher & Oemlerl(ll984l) . Symbols are as in Figures [TOland [TTI We add a new pane l at z — 1.05 for model prediction, where 
we overplot the measured values of fbi ue for the RzCS 052 cluster in Andreon et al.l (120081) . 



redder color to split blue and red galaxies, hence in increasing 

fblue- 

As a test, we repeated our analysis, with a non-evolving, 
mass-dependent definition of blue/red galaxy. For low redshifts 
(z < 0.5), this leads to results marginally different from ours, 
because the colors used in the two methods to split blue and 
red galaxies do not differ significantly. However, for higher 
redshifts, we do observe a notable difference between the two 
methods, because the difference between the two definitions of 
blue/red galaxies will increase with increasing redshift. For in- 
stance, using a non-evolving color threshold leads to larger blue 
fractions at z > 1 (0.2 <, fbiue i5 0.6) for stacked data. Since 
fbiue depends on r/r2oo, a non-evolving color threshold also tilts 
the fbiue profile. Overall, the results using a non-evolving color 
threshold are much harder to interpret, because there is an ad- 
ditional effect that depends on the redshift (and potentially the 
radial) bin(s) under inspection. 



5.2. Remark on using evolved mass definition 

We remark that using the mass in stars at the redshift of obser- 
vation, M zo i, s , as a mass definition does not appreciably change 
our conclusions: the estimated values of fbi m are slightly lower, 
but the trivariate dependence of fbiue on redshift, galaxy mass 
and cluster-centric distance remains similar (the a, /3, y, and 5 
parameters vary by less than lcr). 

Indeed, as illustrated in Figure|8] lines of constant mass at the 
redshift of observation in the color-magnitude diagrams are sim- 
ilar to those used in the study, but for blue colors. For a consid- 
ered mass bin, using M w b s instead of M will remove galaxies in 
the light-gray shaded area, but add those in the dark-gray shaded 
area. Because of the mass function shape, the former will on av- 
erage be slightly more numerous than the latter: the net effect 
will be to slightly decrease the number of blue galaxies, hence 

fblue- 



5.3. Dependence of fbiue on M and z 
The data and their fit point out two results: 

- The average SFH is more complex than exponentially declin- 
ing SFHs, as quantified by the y and f terms. Indeed, if clus- 
ter galaxies had on average exponentially declining SFHs, 
we should not observe any change of the fu ue radial pro- 
file with redshift (at fixed galaxy mass). Low-mass galaxies 
may have a briefly enhanced SFH because of starbust, hence 
temporarily a color bluer than an SFH t - 3.1 Gyr model, 
after which their color turns redder than an SFH t = 3.7 
Gyr model. For our lowest mass bin, this mechanism would 
still be at work at z < 0.37. We also considered the case of 
a later Zform f° r low-mass galaxies. We experimented with 
some possible choices without being able to remove the ob- 
served evolution of fbh, e with redshift, although our tests are 
not exhaustive and we are reaching the limit of the data. 

- There is a differential evolution of fbiue with the galaxy mass, 
i.e. mass quenching is a dynamical process. Galaxies with 
smaller masses evolve later: for instance, at z < 0.37, fbi ue 
for galaxies in our lowest mass bin still evolve significantly, 
which is not the case for fbi„ e for more massive galaxies. 

Furthermore, less massive galaxies seem to evolve on a 
longer timescale than more massive galaxies. In the ~3 Gyr be- 
tween our two lowest redshift bins, fbiue changes by ~0.4 for 
low-mass galaxies. In less time (~2 Gyr), fbiue should change by 
twice as much at z > 2 for high-mass galaxies. 

5.4. Dependence of fbiue on r/r 2 oo 

We now turn to the analysis of the dependence of fbiue with 
cluster-centric distance r/r2oo- For all galaxy stellar mass and 
redshift bins, fbiue decreases with decreasing r/r^ca, i-e. environ- 
mental quenching is effective at all redshifts and stellar masses 
(Figure [TTJ. This effect is conspicuous on lower mass galaxies, 
as fbiue is ver y l° w for massive galaxies, whatever the cluster- 
centric distance (one should observe a significant number of 
clusters at z ~ 1-2 to see this effect on massive galaxies). 



13 



A. Raichoor & S. Andreon: Evolution of galaxies in clusters in the last 10 Gyr 



Previous works clearly establ ished this dependence on r/r2oo for 
mass- selected samples (e.g. lAndreon et al.l l2008t lHaines et alJ 
l2009h . however, our study shows for the first time that this de- 
pendence holds for different galaxy mass bins. In addition, the 
model fitting our data implies that the mass and environmental 
quenchings are fully separable, because no crossed term between 
M and r/r^oo is required. 

Additionally, we found that our data do not require any evo- 
lution of this dependence with redshift: our (simple) model re- 
quires only a crossed term between M and z, none between r/r2oo 
and z- This can be interpreted that either this evolution with z of 
the environmental quenching is not present, or that is of second 
order when compared to the evolution with z of the mass quench- 
ing. We recall that our cluster sample consists of clusters having 
similar temperatur es at all redshi fts, while a self-similarly evolv- 
ing cluster model dKaiser|[l986l) suggests a mild increase with 
redshift at < z < 1 . Nevertheless, current studies tend to show 
that fyiM estimated in the virial radius is independent of cluster 
mass, at lea st in the local Universe (e.g. jDe Propris et al.ll2004l 
lGotoll2005h . 

Finally, we would like to emphasize the importance of the 
backsplash population on the fuue radial profile and its evolu- 
tion with redshift: galaxies as far away as two virial radii may 
have been within the main body of the cluster in the pa st (e.g., 
iBalogh et alJl2000tlGill et a!ll200llMah"aian et al.ll201 ll 

5.5. Conclusion 

Our study extends the relationship between galaxy mass, star 
forma tion rate and environment at z < 1 explored bv lPeng et al.l 
to intermediate-mass cluster environments. As those au- 
thors, we found that the cessation of star formation ("quench- 
ing") due to environment and galaxy mass are separable (we 
needed no crossed term between r/r^oo and M to fit our data), 
that environmental quenching does not change with epoch (we 
needed no crossed term between r/r^oo and z to fit our data), 
and that mass quenching is a dynamical process (we did need a 
crossed term between M and z to fit our data). The close agree- 
ment between the two studies is likely a consequence of the com- 
mon choices adopted for the two analyses: disentangling galaxy 
mass and environment, plus taking into account th e secular aging 
of star s with decreasing redshift, as described by lAndreon et al.l 
However, in addition to paying special attention to the 
statistical aspect of the analysis (which allowed us to provide 
good confidence intervals on our measurements), our study im- 
plemented a finer control on galaxy mass (by considering the 
mass evolved at z — 0) and on environmental estimation (us- 
ing r/r^oo is less subject to biases than a density field based on 
res t-frame fl-band selec tion) . 

iMuzzin et al.l d2012l) . investigating a cluster sample at z ~ 1 
with spectroscopic data, made a similar analysis, studying the 
role of galaxy mass and environment (but not secular evolution) 
on star formation. The availability of spectroscopic data allows 
estimating the star formation rate, and thus a thorough analysis 
of this problem. Regarding the fraction of star-forming galax- 
ies, these authors also concluded that mass and environmental 
quenching are separable. Our work extends this study by includ- 
ing the redshift dependence in the analysis, by a better character- 
ization (f2oo) of cluster properties, and by using a cluster sample 
that is not selected by a galaxy property under study (color). 

Our study extends the downsizing-like scenario to cluster 
environment and all galaxies, which was already established 
for galaxi es with specific p roperties in the field (e.g., emis- 
sion line: ICowie et alJll996t spheroidal: iTreu et al 1 I2001 pas- 



sive: iPeng et ai1l2010l) . According to this scenario, the proper- 
ties of the most massive galaxies are established in the very 
early Universe (z » 1), while less massive galaxies continue 
to evolve at redshifts < z < 1. We stress that our down- 
sizing-like scenario concerns all ga laxies, regardless of their 
morphology or color. The r esults of lAndreon et al.l (120081) and 
iRaichoor & Andreonl (1201 2l) taken as a whole point to a similar 
conclusion, though for smaller cluster samples. 

Our work is the first to study the dependence of fbi ue in clus- 
ters on galaxy mass, redshift, and cluster-centric distance at the 
same time, with a cluster sample whose properties are well- 
controlled. In particular, it clearly demonstrates the need to use 
galaxy mass as a parameter to better understand the behavior of 
fbiue ■ This approach is now possible with the wealth of available 
data. 
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Appendix A: r 2 oo estimation 



measurements dPacaud et all l2007t 
is determined using the mass- 



For clusters with Tx 
iHudson et all 1201 Oh . r ?m 

temperature relation of lFinoguenov et all (1200 ll) . which has been 
shown to give ro bust results for clusters with Tx < 4 keV 
(IWillis et al.ll2005l) : 



r 500 (Mpc) 



0.391 • r°- 63 
H(z)P0 ' 



(A.l) 



where H(z) = Ho -\j£l m (l + z) 3 + Qa describes the redshift 
evolution of the Hubble parameter in the assumed cosmological 
model. Then we us e r2oo = '"500/O.66I corresponding to a 
iNavarro et al. I (11997b profile with a halo concentration parame- 
ter c = 5. 

For clusters with no Tx measurements (XLSSC007, 
XLSSC014, and XLSSC016), we estimate r 2 po fro m the clus- 
ter velocity dispersion cr v , from IWillis et all d2005l) . following 
the formula 



r 20 o (Mpc) 



V3 • <r v [kms '] 
10 ■ H(z) 



xO.85. 



(A.2) 



This equation is similar to the classical equatio n for a spher- 
ical collapse model (e.g., ICarlberg et al.l 1 19971) . but for the 
multiplicative coefficient 0.85. This coefficient comes from a 
calibration we perform with a compilati on of clusters havin g Tx 
and cr v measurements (see for instance [Proctor et al.ll201 ll for a 
similar systematic difference). Eq. (lA.2b ensures that there is no 
bias between our r2oo estimated with Tx or <x v (see Figure lATTT i. 
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Fig. A.l. Comparison of rooo estimated with Eqs. flA.ll i & 1A.21 : 
the test sample is made of 10 clusters from the XMMLSS sur- 



vey (IWillis et al.ll2005HAndreon et al.ll2006l:lPacaud et al.l2007 | ) 
28 clusters from the CAIRNS su rvey (iRines & DiaferiolhoOol 
and 15 clusters from the CNOC 1 dMuzzin et al.ll2007l) and spans 
0.003 < z S pec < 0.84, 0.5 < T x (keV) < 10.3 and 232 < cr v (km 
s~') < 1354. There is no bias between the r2oo estimated from 
Tx and from <x,,. The horizontal (resp. vertical) arrows indicate 
the riot) estimated from Tx (resp. cr v ) for our cluster sample. 




20 

r-band magnitude 

Fig. B.l. zo and y parameters used for building r-band prior with 
Eq.(EB. 



as in lBrammer et al.l (120081) . At 13.5 < r and 18 < r < 20, we 
visually extrapolated zo and y (cf. Figure IB. It . because either 
there are not enough galaxies or those present in the SDSS are a 
biased sample (e.g., LRGs, QSOs, etc). 



Appendix B: r-band prior for photometric redshifts 

Because Eazy has been developed for high-z studies, the r-band 
prior does not include magnitudes brighter than 20. In our SDSS 
sample, we deal with galaxies as bright as r ^ 12. We built the 
prior probabilities p(z\mo) for 13.5<r<18as follows. For each 
interval of 0.5 mag, we retrieved ~ 10,000 spectroscopic galaxies 
from the SpecPhoto view from the SDSS DR8 database and 
fitted the redshift distribution with a function 
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p(z\m ) oc z y ■ exp [-(z/zo) 7 ], 



(B.l) 
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